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ABSTRACT 

Using direct simulations, weakly nonlinear theory and nonlinear mean-field theory, it is 
shown that the quenched velocity field of a saturated nonlinear dynamo can itself act 
as a kinematic dynamo. The flow is driven by a forcing function that would produce a 
Roberts flow in the absence of a magnetic field. This result confirms an analogous find- 
ing by F. Cattaneo & S. M. Tobias (arXiv:0809.1801) for the more complicated case of 
turbulent convection, suggesting that this may be a common property of nonlinear dy- 
namos; see also the talk given also online at the Kavli Institute for Theoretical Physics 
(http://online.kitp.ucsb.edU/online/dynamoj:08/cattaneo). It is argued that this property can 
be used to test nonlinear mean-field dynamo theories. 
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1 INTRODUCTION 

The magnetic fields of many astrophysical bodies displays order 
on scales large compared with the scale of the turbulent fluid mo- 
tions that are believed to generate these fields via dynamo action. A 
leading theory for these types of dynamos is m ean-field electrody- 
namics <Moffadll973 : lKrause & Radlerlll980h . which predicts the 
evolution of suitably averaged mean magnetic fields. Central to this 
theory is the mean electromotive force based on the fluctuations of 
velocity and magnetic fields. This mean electromotive force is then 
expressed in terms of the mean magnetic field and its first derivative 
with coefficients ay and 77yfe. The former represents the a effect 
and the latter the turbulent magnetic diffusivity. 

Under certain restrictions the coefficients ay and rjijk can 
be calculated using for example the first order smoothing approx- 
imation, which means that nonlinearities in the evolution equa- 
tions for the fluctuations are neglected. Whilst this is a valid ap- 
proach for small magnetic Reynolds numbers or short correla- 
tion times, it is not well justified in the astrophysically interest- 
ing case when the magnetic Reynolds number is large and the cor- 
relation time comparable with the turnover time. However, in re- 
cent years it has become possible t o calculate aq and r / qk us - 
ing the so-called test-field method dSchrinner et alj|2005l . 120071) . 
For the purpose of this paper we can consider this method essen- 
tially as a "black box" whose input is the velocity field and its out- 
put are the coefficients ay and rjijk- This method has been suc- 
cessfully applied to the kinematic case of weak magnetic fields 
in the presence of homogeneous turbulence either witho ut shear 



dSur et al.ll2008l; IBrandenburg. Radler & Schrinnerl l2008h or with 
shear dBrandenburg||2005h IBrandenburg et all2008al^ 

More recently, this method has also been applied to the non- 
linear case where the velocity field is modified by the Lorentz force 
associa ted with the dynamo-generated field ( IBrandenburg et all 
2008b) In that case the test-field method consists still of the same 



black box, whose input is only the velocity field, but now this ve- 
locity field is based on a solution of the full hydromagnetic equa- 
tions comprising the continuity, momentum, and induction equa- 
tions. We emphasize that the magnetic field is quite independent of 
the fields that appear in the test-field method inside the black box. 

Our present work is stimulated by an intere sting and relevant 
numer ical experiment performed recently by ICattaneo & Tobiaa 
(2008). They considered a solution of the full hydromagnetic equa- 
tions where the magnetic field is generated by turbulent convective 
dynamo action and has saturated at a statistically steady value. They 
then used this velocity field and subjected it to an independent in- 
duction equation, which is equal to the original induction equation 
except that the magnetic field B is now replaced by a passive vec- 
tor field B, which does not react back on the momentum equation. 
Surprisingly, they found that |J3| grows exponentially, even though 
the velocity field is already quenched by the original magnetic field. 

One might have expected that, because the velocity is modi- 
fied such that it produces a statistically steady solution to the origi- 
nal induction equation, B should decay or also display statistically 
steady behavior. The argument sounds particularly convincing for 
time-independent flows because, if a growing B were to exist, one 
would expect this alternative field to grow and replace the initial 
field. This view is supported by recent simulations in which the 
flow field from a geodynamo simulation in a spherical shell was 
used as velocity field in kinematic dy namo computations, and no 
growing B was found dTilgnerl l2008). However, it turns out that 
this reasoning is not correct in general. One finds counterexamples 
even within the confines of mean field MHD using analytical tools. 
The existence of a growing B thus is not tied to chaotic flows or 
fluctuating small-scale dynamos. 

This finding of ICattaneo & Tobias! (2008) is interesting in 
view of the applicability of the test-field method to the nonlinear 
case. Of course, the equations used in the test-field method are dif- 



© 2008 RAS 



2 A. Tilgner & A. Brandenburg 



the context of a simpler example. However, let us begin by repeat- 
ing the numerical experiment of ICattaneo & Tobias! J2OO8) using 



ferent from the original induction equation. (The equations used in 
the test-field method include an inhomogeneous term and the mean 
electromotive force is subtracted, but they are otherwise similar to 
the original induction equation.) Given the seemingly unphysical 
behaviour of the induction equation in the presence of a vector field 
different from the actual magnetic field, it would be tempting to ar- 
gue that one should choose test fie lds whose shape is rathe r close 
to that of the actual magnetic field iCattaneo & H ughes 200^). On 
the other hand, the a%j and rjijk tensors should give the correct re- 
sponse to all possible fields, not just the B field that grew out of a 
particular initial condition, but also the passive B field that obeys a 
separate induction equation. It is therefore important to choose a set 
of test fields that are orthogonal to each other, even if none of the 
fields are solutions of the induction equation. One goal of this paper 
is to show that the ctij and rjijk tensors obtained in this way pro- 
vide not only interesting diagnostics of the flow, but they are a lso 
able to explain the surprising result of lCattaneo & Tobia s N2008h in 

» begin by repea 
3lasN2008l) usir 

the simpler case of a Roberts flow. Next, we consider a weakly 
nonlinear analysis of this problem and turn then to its mean-field 
description. 



2 THE MODEL 

In order to examine the possibility of a growing passive vector field, 
we first considered the case of a driven ABC flow. Such a flow is 
non-integrable and has chaotic streamlines. Growing passive vec- 
tor fields were found. To simplify matters even further, we consider 
now the case of a Roberts flow, which is integrable, has non-chaotic 
streamlines, and the dynamo can only be a slow one, i.e. the growth 
rate goes to zero in the limit of large magnetic Reynolds number. 
This is however not an issue here, because we will only be consid- 
ering finite values of the magnetic Reynolds number. 

In the following we consider both incompressible and isother- 
mal cases. The governing equations for any externally driven ve- 
locity field (turbulence, ABC flow, or Roberts flow) are then given 
by ' 

^- = -U-VU-VH+-JxB + f + F visc , (1) 
at p 

^ = Vx([/xB)+7jV 2 B, (2) 

where U is the velocity, B is the magnetic field, p is the density, 
H is the specific enthalpy, J — V x B/po is the current density, 
po is the vacuum permeability, / is the forcing function, -Fvisc is 
the viscous force per unit mass, and 77 = const is the magnetic 
diffusivity. In the incompressible case, V ■ U = 0, we have H = 
p/p, where p is the pressure and p = const. The viscous force is 
then given by .Fvisc = vV 2 U. In the isothermal case, the density 
obeys the usual continuity equation 

|=-V-(pt0, (3) 

but now p — c 2 p, where c s is the isothermal sound speed. In that 
case H = c 2 In p and the viscous force is given by 



F visc = vV 2 U + |^VV ■ U + 2uSV ln(pu), 



(4) 



where Sy = \(Ui,j + t/j',») — 3&jV • U is the traceless rate of 
strain matrix. 

In order to compute the evolution of an additional passive vec- 
tor field B we also solve the equation 



^=Vx([/xB)+ V V 2 B. (5) 

at 

In the case of the Roberts flow we use the forcing function 

/ = vk 2 U Rob, (6) 

where 

Uroo = k t ipz -zxVip (7) 
with 

ip = (uo/ko) cos kox cos koy (8) 

and kf = y2ko. We consider a domain of size L x x L y x L z . In all 
cases we consider L x = L y = L z = 2n/ko. Our model is charac- 
terized by the choice of fluid and magnetic Reynolds numbers that 
are here based on the inverse wavenumber ko, 

Re = u Q /vk , R m = u /r)k . (9) 



3 NUMERICAL EXPERIMENTS 

We solve equations {TJ-((5} for the isothermal case using the PEN- 
CIL CodfJj, which is a high-order public domain code (sixth 
order in space and third order in time) for solving partial dif- 
ferential equations. Equation 101 is solved using the test-field 
module with the input parameters lignore_uxbtestm=T, and 
itestfield='B=0', which means that the inhomogeneous 
term of the test-field equation is set to zero and the subtraction 
of the mean electromotive force has been disabled. In this way 
we solve equation $5$, instead of the original test-field equation. 
We focus on the case of small fluid Reynolds number, Re = 0.5. 
The initial conditions for B and B are Beltrami fields, (cos(fco£ + 
ifi), sin(fcoz + <p), 0), with an arbitrarily chosen phase ip — 0.2, but 
for B we ip = 0. In Fig. Q] we plot the evolution of the rms val- 
ues of B and B for a weakly supercritical case with R m = 6.25. 
(In our case with Re = 0.5 the critical value for dynamo action is 
Rm ~ 5.77; for Re — > the critical value would be i? m « 5.5.) 
Both B and B grow at first exponentially at the same rate. How- 
ever, when B reaches saturation, the growth of B slows down tem- 
porarily, but then resumes to nea r ly its o riginal value. This confirms 
the result of ICattaneo & Tobias! J2008b for the much simpler case 
of a Roberts flow. 

In Fig. [2] we compare horizontal cross-sections of the two 
fields. Note that the two are phase shifted in the z direction by a 
quarter wavelength. The short interval in Fig. Q] during which the 
growth of B slowed down temporarily is therefore related to the 
fact that the solution needed to "adjust" to this particular form. 



4 WEAKLY NONLINEAR THEORY 

A weakly nonlinear an alysis of the Roberts flow is presented in 
iTilgner & Bussd u200lh . Two nonlinear terms enter the full dynamo 
problem. In order to make the calculation analytically tractable, it 
is assumed that the fluid has infinite magnetic Prandtl number so 
that the inertial terms (and hence the advection term) drop from the 
Navier-Stokes equation. The second nonlinear term, the Lorentz 
force, is assumed to be small compared with the driving force / 
and is treated perturbatively. The linear induction equation is trans- 
formed into a mean field equation assuming separation of length 



1 jhttp: / /www, nordita . org/ software /pencil- code /| 
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Figure 1. Evolution of the rms values of B and B for ij m = 6.25. The 
growth rate of B is 0.016 itofco m me kinematic phase and 0.014 Ufjfco in 
the nonlinear phase. 



scales and small magnetic Reynolds numbers. One can under these 
assumptions compute the modifications of the velocity field Unob 
due to the presence of a mean field B, which we define here as 



B(z,t)= / B{x,y,z,t)dxdy/L x L y . (10) 

Jo Jo 

The magnetically modified velocity field then becomes 

(l-7)^Rob 



U 



with 



U 



Uo 



B,. + B 



sin kox cos koy N 
— cos kox sin koy 
, \f2 sin kox sin koy , 

/-ft2 -=;2 , 



U 



(ID 



(12) 



and 7 = (£>^ + By)/(2r]vkj ppo)- The original flow t/R b is 
reduced in magnitude and another 2D periodic flow component is 
added. The mean-field induction equation with flow U given by 
i ll It . written for a passive vector field B, is: 



dB 

dt 

with 



V x A 




- V x CB X B 




r]V 2 B 



(13) 



A = R m vo and C - 



B'r, 



1 



where 

5 _ v 

T)kf 



p - V _ 



and «o = |wq(1-7) 



(14) 



(15) 



ha ve been introduced. Thi s equation corresponds to equation (10) 
of lTilgner&Bussd l l200ll) . Apart from a change of notation, the 
distinction between the passive vector field B and the field distort- 
in g the Roberts flow, B, h as been made. In addition, equation (10) 
of iTilgner & Bussd feOOll) was intended as a model of the Karls- 
ruhe dynamo and used uo as control parameter, whereas here, we 
consider / as given. For this reason, the two equations are identical 
only to first order in 7. 

Equation U3b reduces to the usual dynamo problem for B = 
B and leads to the kinematic dynamo problem if it is further- 
more linearized in B, which corresponds to dropping the third 
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Figure 2. Grey-scale representations of horizontal cross-sections of B and 
B for ij m = 6.25. Here, kozo ~ —3.04. Both fields are scaled symmetri- 
cally around zero (grey shade) with dark shades indicating negative values 
and light shades indicating positive vales, as indicated on the grayscale bar. 



term in Eq. d!3t . For uq > 2rj 2 k{k, the equation has growing so- 
lutions of the form (cos kz, sin kz, 0). Weakly nonlinear analysis 
determines the amplitude Bo of a saturated solution of the form 
B — Bo(coskz,sinkz,0) by inserting this ansatz for B and 
B = B into equation d 1 3 b and by projecting equation d 1 3b onto 
(cos kz, sin kz, 0) and integrating spatially over a periodicity cell 
( ITilgner & Bussell200lh . 

We proceed similarly to find growing solutions of equa- 
tion dl3t . Assume B — £?o(cos kz, sin kz, 0). Obvious candidates 
for passive vector fields growing at rate p are B oc e p (cos(kz + 
95), sin(fcz + ip), 0). Since the third term in equation d 1 3 b is a per- 
turbation, a growing solution must be of a form such that the other 
terms maximize the growth rate. These other terms are identical to 
the kinematic problem, so that B must have the same general form 
as the kinematic dynamo field except for the phase shift ip which 
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measures the phase angle between the saturated field B and the 
passive vector field B. The velocity field Unoh is independent of z, 
so that any solution of the kinematic dynamo problem remains a so- 
lution after translation along z. However, neither the Lorentz force 
due to B nor the flow modified by that Lorentz force are indepen- 
dent of z, so that the phase angle tp matters for B. The above ansatz 
for B will not be an exact solution of equation d 1 3 fc because of the 
z-dependence of B x B y , but it represents the leading Fourier com- 
ponent. In order to determine the optimal tp, we insert this ansatz 
into equation J 1 3 1 . project onto (cos(fcz + tp), sin(fcz + </?), 0) and 
integrate over a periodicity cell to find 



p — —rjk + Ak 



CBgk— cos 2ip. 



(16) 



The saturation amplitude Bo is determined from this equation by 
setting p = and tp = 0. For any given Bo, the maximum of p 
is obtained for tp = n/2. The fastest growing mean passive vector 
field is thus expected to have the same form as the mean dynamo 
field except for a translation by a quarter wavelength along the z- 
axis. This is in agreement with the simulation results of Sect. [3] 

The weakly nonlinear analysis in summary delineates the fol- 
lowing physical picture: As detailed in Tilgner & Busse ( 200lj), the 
saturated dynamo field modifies the flow in two different ways. 
Firstly, it reduces the amplitude of JJ-^oh by the factor 1 — 7, and 
secondly, it introduces a new set of vortices which lead to the third 
term in equation d!3t . The reduction of the amplitude of the Roberts 
flow affects all magnetic mean fields with a spatial dependence in 
(cos(fcz + tp), sin(A;2; + tp), 0) in the same manner, independently 
of tp. The additional vortices, however, have a quenching effect on 
the field that created them (e.g. tp = 0) but are amplifying for a 
field shifted by tp — tt/2 with respect to the saturated field. 

We were able to find a simple growing passive vector field 
thanks to the periodic boundary conditions in z. The same construc- 
tion is impossible for vacuum boundaries at z = and z = 2-n/ko- 
Numerical simulations of the Roberts dynamo with vacuum bound- 
aries, not reported in detail here, have revealed that growing B ex- 
ist in this geometry nonetheless, but they bear a more complicated 
relation with B than a simple translation. At present, the flow of 
the co nvection driven dynamo in a spherical shell used in [Tilgner 
(2008) seems to be the only known example of a dynamo which 
does not allow for growing B. 



5 NONLINEAR MEAN-FIELD THEORY 

In mean-field dynamo theory for a flow such as the Roberts flow 
one solves an equation for the horizontally averaged mean field, as 
defined in equation dlO| >. The mean electromotive force is defined 
as £ = u x b, where u — U — U and b = B B are the 
fluctuating components of magnetic and velocity fields. The mean 
electromotive force can be expressed in terms of the mean fields as 



rjij J j 



(17) 



where we have used the fact that for mean fields that depend only 
on one spatial coordinate one can express all first derivatives of the 
components of B in terms of those of J alone. 

The tensorial forms of otij and rjij are ignored in many mean- 
field dynamo applications, but here their tensorial forms turn out to 
be of crucial importance. Of course, there is always the anisotropy 
with respect to the z direction, but this is unimportant in our one- 
dimensional mean field problem, because of solenoidality of B and 
J and suitable initial conditions on B such that B z = J z — 0. 



However, the dynamo-generated magnetic fields will introduce an 
anisotropy in the x and y directions. If B is the only vector giving 
a preferred direction to the system, the ay and rjij tensors must be 
of the form 



an = ai(B)8ij + a 2 (B)B t Bj + 
rjij = vi(B)Sij + rj 2 (S)BiBj + 



(18) 
(19) 



where B — B/\B\ is the unit vector of the dynamo-generated 
mean magnetic field, and dots indicate the presence of terms related 
to the anisotropy in the z direction inherent to the Roberts flow. As 
indicated above, equation d 1 8b is correct without these terms only 
in the (a;, y) plane. However, the terms represented by the dots do 
not enter the considerations below because we are only interested 
in fields with B z — 0. 

In order to predict the evolution of B in the saturated state, 
we need to know the effect on ctij (and in principle also on rjij, 
but 772 is small; see below). Thus, we now need to know B. The 
mean magnetic field generated by the Roberts flow is a force-free 
Beltrami field of the form 



B — (cos koz, sin koz, 0), 



BiBj = 



cos koz 
cos koz sin koz 




cos koz sin koz s 
sin 2 koz 



(20) 



(21) 







0, 



The coefficients ai, 012, and r/2 have previ ously been deter- 
mined for the case of homogeneous turbulence dBrandenburg et all 
2008b) and it turned out at oti and a.2 have opposite sign, and that 
r)2 is negligible. This is also true in the present case, for which we 
have determined cti/uo = —0.266, 012 juo = +0.022, rjiko/uo = 
0.082, and rjzko/uo = 0.002. The microscopic value of 77 is 0.160, 
so the steady state condition, ai + a.2 + {fji + r\2 + r])ko = 0, 
is obeyed[f| In the kinematic regime we have ai/ito = —0.254, 
rjiko/uo = 0.076, with 012 = r\2 — 0, resulting in a positive 
growth rate of 0.018 uoko- Thus, even though ai increases in this 
case, the sum ai + a?2 is being quenched. This, together with the 
increase of 771 + 772, leads to saturation of B. 

Returning now to the mean-field problem for B, this too will 
be governed by the same ctij and rjij tensors, but now the tensor 
BiBj is fixed and independent of B. The solution for B will be one 
that maximizes the growth, so it must experience minimal quench- 
ing. Such a solution is given by that eigenvector of BiBj that min- 
imizes the quenching of B. In the case of our Beltrami field (HQ, 
the minimizing eigenvector is given by 



B — (sin koz, — cos koz, 0), 



(22) 



0. This is indeed the same result that 
we found both numerically and using weakly nonlinear theory. The 
growth rate of B is then expected to be |«i|fco — + Tj)^o = 
0.024 uoko, which is indeed positive, but it is somewhat bigger than 
the one seen in Fig.Q] 

Let us emphasize once more that by determining the full a%- 
and rjij tensors in the nonlinear case, we have been able to predict 
the behavior of the passive vector field as well. This adds to the 



2 We note that the sign of a\ is opposite to the sign of the kinetic helicity, 
but since the Roberts flow has positive helicity, Qi must be negative, which 
is indeed the case. 
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credence of the test-field method in the nonlinear case, and con- 
firms that the test-fields can well be very different from the actual 
solution. 

The considerations above suggest that solutions to the passive 
vector equation equation ((5} can be used to provide an independent 
test of proposed forms of a quenching. Isotropic formulations of a 
quenching would not reproduce the growth of a passive vector field, 
and so such quenching expressions can be ruled out, even though 
the resulting electromotive force for B would be the same. We sug- 
gest therefore that the eigenvalues and eigenvectors of equation l[5j 
with a velocity field from a saturated dynamo can be used to charac- 
terize the quenching of dynamo parameters (a effect and turbulent 
diffusivity) and thereby to test proposed forms of a quenching. 



6 CONCLUSIONS 

The most fundamental question of dynamo theory beyond kine- 
matic dynamo theory is "how do magnetic fields saturate?". In the 
simplest picture, the velocity field reorganizes in response to the 
Lorentz force such that all magnetic fields decay except one which 
has zero growth rate and which is the one we observe. This pic- 
ture is already questionable for chaotic dynamos. In a chaotic sys- 
tem, nearby initial conditions lead to exponentially separating time 
evolutions. If one takes a magnetic field B with (on time average) 
zero growth rate which is the saturated solution of a chaotic dy- 
namo, and solves the kinematic dynamo problem for a passive vec- 
tor field B with initial conditions different from B, one is prepared 
to find growing B. Examples for growing B in chaotic dynamos 
have been given bv lCattaneo & Tobias! |2003) . 

For a time-independent saturated dynamo, on the other hand, 
the simple picture seems to be adequate at first. However, we 
have shown in this paper that growing B also exist in the time- 
independent Roberts dynamo. The origin of the growing B can in 
this case be understood with the help of weakly nonlinear theory. 
The growing B has the same shape as the saturated dynamo field 
but is translated in space. 

What was wrong with the naive intuition invoked above? It 
was based on a stability argument (the equilibrated magnetic field 
should be replaced by B if there is a growing B). However, the 
linear stability problem for a solution of the full dynamo equations 
is different from the kinematic dynamo problem for B, because in 
the latter, the velocity field is fixed. Both problems are closely re- 
lated eigenvalue problems, but standard mathematical theorems do 
not provide us with a relation between the spectra of both prob- 
lems. The numerical computation above gives an example of a sta- 
ble Roberts dynamo, showing that the linear stability problem for 
the solution found there has only negative eigenvalues. Solving the 
same eigenvalue problem with velocity fixed, which is the kine- 
matic problem for B, can very well lead to positive eigenvalues, 
and indeed, it does. The dynamo is therefore only stable because 
the velocity field is able to adjust to perturbations in the magnetic 
field. The magnetic field on its own is unstable. 
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